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Methods for Correcting Microwave Scattering and 
Emission Measurements for Atmospheric Effects 

by Mark Komen 


ABSTRACT 

Algorithms have been developed to permit correction of scattering coefficient 
and brightness temperature for the Skylab S"I93 Radscat for the effects of cloud 
attenuation. These algorithms depend upon a measurement of the vertically polarized 
excess brightness temperature at 50° incidence angle. This excess temperature is 
converted to an equivalent 50° attenuation, which may then be used to estimate the 
horizontally polarized excess brightness temperature and reduced scattering coefficient 
at 50°. ‘^''r angles other than 50° the correction also requires use of the variation 
of emissivity with salinity and water t- inperature. When theoretical values are used 
for this, the effect may be translated to an equivalent 50° case and the attenuation 
estimated. 

Use of 50° estimates should be quite satisfactory for the in-track modes of the 
S-193, but the more complex (and less reliable) methods based on using theory to obtain 
an equivalent 50° number must be used in the cross-track modes because the same spot 
is not viewed at both 50° and another angle in these modes. 


1.0 INTRODUCTION 


The removal of atmospheric effects from radiometer and scatterometer measure” 
ments has been the focus of this part of the technical study on the RADSCAT. Two 
main areas, correction of scattering coefficient for atmospheric effects and brightness 
temperature for atmospheric, angle, and surface effects are discussed in this memo” 
randum. 

It has been shown (Hollinger, 1971) that vertically polarized brightness temp- 
erature at 50® incidence angle is insensitive to surface roughness and correspondingly 
to windspeed while still being sensitive to changes in the atmosphere. Brightness 
temperature at the surface is related to apparent temperature at the satellite altitude, 
the difference bein ed as excess temperature. If the excess temperature at 50° 

incidence angle is known, the theory covering the 50° incident angle case can be ex- 
tended to all incidence angles. 

Also to be considered is translating the brightness temperature to a reference 
brightness temperature at some standard water temperature and standard incidence angle 
so that the radiometric measurements after the removal of the temperature and incident 
angle effects reflect only atmospheric and roughness effects. The final result of correct- 
ing for atmospheric effects is presented in the form of two user-oriented subroutines, one 
which removes the atmospheric effects based on the 50° incident angle brightness 
temperature insensitivities, and the other which performs the reference temperature/ 
incidence angle translotion. 

The purpose of this memorandum is twofold, the first being to remove the effects 
of the atmosphere from the radiometer and scatterometer measurements. 

it has been shown (Hollinger, 1971) that brightness temperature at 50° incident 
angle Is insensitive to surface roughness and correspondingly to windspeed while still 
being sensitive to changes in the atmosphere. Brightness temperature at the surface is 
related to apparent temperature at the satellite altitude, the difference beirig defined as 
excess temperature. F the excess temperature at 50° incident angle is known, the 
theory covering the 50° incident angle case can be extended to all incident angles. 

Also to be considered Is translating the brightness temperature to a reference 
brightness temperature at some standard water temperature and standard incident angle 
so that the radiometric measurements reflect the removal of the temperature and incident 
angle effects. The final result of correcting for atmospheric effects is presented in the 


form of two user-orlented subroutines, one which removes the atmospheric effects based 
on the 50° incident anjie brightness temperature insensitivities, and the other which 
performs the reference temperature/incid'^nt angle translation. 

2.0 Removing the Effects of the Atmosphere~Derivotion of Excess Temperature and 
Attenuatio n 

The estimation of the atmospheric effect begins with the calculation of vertical 

excess temperature near 50° incidence angle. This excess temperature is calculated by 

taking the difference in the apparent temperature measured by Skylab and the surface 

brightness temperature as calculated from Fresnel theory. Although Fresnel theory is 

only valid for smooth surface conditions at 50° incidence angle and vertical polarization, 

the result is independent of windspeed. Given the incident angle, salinity, frequency, 

and the surface water temperature, the reflection coefficients for an electromagnetic 

wave incident on a smooth sea (zero wind case) can be calculated. The polarized 

reflection coefficient R (6,Tw), where 0 Is the incident angle, Tw is the water temperature, 

P 

and p represents the polarization, is expressed as: 




CC.S0 

6v-(Tw) 


^'0 


where (Tw) is the complex relative dielectric constant of sea water at temperature Tw. 
Emissivity ep (6, Tw) is related to the Fresnel coefficient by 


^ I - I l?p 








0 ) 


The vertical brightness temperature is expressed as: 

Tb, = Tw)% 

Where is the water temperature at the surface. By definition, the excess temperature 
is the difference between the apparent temperature and the brightness temperature. 

where T is the radiometric temperature at the satellite altitude in the antenna 
app 

pointing direction and Tg is the surface brightness temperature. 

An empirical relation between vertical excess temperature and attenuation at 50° 
incidence angle was derived from soundings taken by a NIMBUS satellite. Further 
dialysis by the ATRAD* atmospheric radiation package yielded a graph of attenuation 
vs. vertical excess temperature (Figure 1). Regression analysis shows the two quantities 
related by the equation 

Z V/ 3 


where X represents the vertical excesf temperature in degrees Kelvin andr<(50) the 
attenuation In decibels. 

ATRAD also yielded a relation between horizontol and vertical excess temperatures 
at 50° incidence angle which (Figure 2) by using a regression analysis, is described by 


~ |.53l \/ f \/ — o.oooo? \/ 

where V is the vertical excess temperature and T^^j^ is the horizontal excess temperature 
at 50°incident or !e, all in degrees Kelvin. 

Since the a’ ation by the atmosphere is now known, the scatterometer measure"* 
ment is correcte 
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Figure 1. Attenuation vs. vertical excess temperature (50® incidence angle). 



Figure 2, Horizontal excess temperature v$, vertical excess temperature 
(50*^ Incidence angle). 




where cr is ^he scatterometer cross section,'^ the attenuation, ando-^ ^he corrected 

c 

scatierometer measurement, all in decibels. 

At this point, all the 50® incidence angle data can be corrected for atmospheric 
effects. Knowing the horizontal excess temperature, the horizontal brightness temp- 
erature at 50® incidence angle can be obtained from equation (2). The horizontal 
emissivity at 50® can be obtained by solving equation (1) for for horizontal 

polarization. 

The next step Is to calculate certain parameters to be used at other incidence 
ongles. 

Attenuation and transmit; jnce are related by the expression 


r(B) = 10 
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wherer(O) is the transmittance at some incidence angle and 'x'' is the attenuation in 
decibels. Transmittance is also a function of incidence an^jle (0) and atmospheric 


opacity (r^). 


-T„sec© 


r(e)= 


Solving for IT f which is independent of incidence angle: 


T = 






It is important to examine the actual compon its of the apparent temperature. Figure 
3 shows an electromagnetic wave incident on some surface having some brightness 
temperature T^, Tg represents the amount of energy the antenna will see transmitted 
through the atmosphere. T^p is the atmospheric contribution T^j^^ (looking up) plus that 
part of the cosmic background T^ (2.7® Kelvin) transmitted through the atmosphere. 
There will be some reflectivity at the surface, so represents that part of T^p 

which is reflected and transmitted. Finally, there is some contribution frofri 

just the atmosphere itself. Summing these components, it is seen that 


Tapp, ^ ^ pT'.p^ 
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SubsfitutinguT^ for Tg, for ^ and (T^j^y 
equation (6) becomes 



+ ri^) for Typ, where 


showing the angular dependence, or 



Solvip^. for (50°) 


■' f(so”)[i-Sv(so‘r«M + i 


or the mean atmospheric temperature is a weighted average of the physical 
temperature of the atmosphere as it variei with altitude and is related to the atmospheric 
temperature and the transmittance (f) by the followir.g expression (see Appendix 

A"2 for o detailed development of 
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is independent of incident angle for a given scan. 

As was shown in equation (4) and (5), given the opacity at 50° and the incident angle, 
the transmittance at any angle can be calculated. Hence, the attenuation at any angle 
is given by 


o((e) = - 


Now, the scatterometer can be corrected as per equation (3). 

Solving for (P)in equation (9) for any incidence angle, and with 

constant for o given scan, 

T„„W*t4,-r(e)] 

Solving for the emissivity at any angle (6, in equation (7) knowir^ the 
transmittance, apparent temperature, water temperature, and the atmospheric tern* 
perature at any incidence angle. 

Since equations for the brightness temperature and excess temperature have been 
defined in equations (1) & (2), the data for incidence angles other than 50® can be 
corrected for atmospheric effects. 

3.0 Removing the Effects of Water Temperature ond lncid '» nt Angle Variations 

To illustrate variation of brightness temperature with windspeed, and, indeed, 
to provide a suitable reference for attet ation determination, effects of other 
parameters must be accounted for. Skylab data were taken at many water tempera- 
tures and incidence ungles; variations in these parameters affect the measured 
radiometric brightness temperature. The incidence angles varied by two or three 
degrees at each command angle while the sea water temperatures varied by 10 or 15 
degrees Kelvin. A nominal seo water temperature of 290®K and nominal incidence 
angles (50.0°, 40.0^ \ 30.0®, 15.0®, and 0.0®) for each command angle were chosen 
and the brightness temperature which would have been measured at these nominal 
conditions was e:,timated. The estimated brightness temperature for these conditions 
is refetTed to as the reference brightness temperature, 

Emissivities for rough surfaces are larger than those for smooth surfaces, but 
Fung and Claassen ( 1] showed that for a given windspeed and frequency, the emis- 
sivities for smooth surfaces and rough surfaces differed by a value which was indepen- 
dent of water temperature. Therefore, the change in emissivity in translating from the 
surface water temperature to the standard water tempercture for a smooth sea will be 


identical with that for a roughened sea. 

Fresnel th. -^ry was used to obtain the pol irized emissivities for a smoorh sea 

os was done previously, fiven frequency, salinity, incidence angle, ond surface 

water temperature (denoted by the point c^, on the W^=0 curve in Figure 4). 

The polarized emissivities are then calculated at the given salinity ond frequency 

o 

but at the standord water temperature of 290 Kelvin (denoted by the point (e^ Tj) 
on the W^=0 curve in Figure 4), The polarized emissivities corrected for atmospheric 
effects (section I) are reorcr-ented by the point (^ 2 ^ ^be curve in Figure 4, 

where W| is the windspetd at the surface cell. Translating to point T^) is done 

by noting 


' w 

Figt 5 shows the emic 'vity vorying with incidence ongle for a given wind" 

speed. It is not known at this time exactly what the incidence angle dependi nee 

is for all windspeeds due to a lack of experimental data. At some points, ine 

variation in windspeed caused by variations in incidence angle are significant and 

therefore the W curv. was used to estimate some first order correction for the data, 
o 

This correction for incidence angl variation Is obtained by translating along the 
curve and assuming the! this is comparable to n-anslating along the curve. 

The total change In emissivlty is then, 

— /I <^. J . +<^16 

Hence, 
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Figure 4, Emissivity vs. water leitiperaK! -t o wind rougher'**d (W|) 
or n dead calm (Wq) sea at o ant incidence angle. 





Knowing the emissivity at the standard water temperature and the standa ■! 
incidence angle, the reference brightness temperature is calculated and the radio- 
meter is corrected for variations in water temperature and incidence !e. 

where Tgi^^p is the reference brightness temperature ond (©^, the 

polarized emissivity at the standard incidence angle / and the standard water 

s 

temperature 
4,0 Algorithms 

This section outlines in flow chart form the algorithms for computing the 
attenuation and emissivities at 50° incident angle (subroutine ATM50), those at any 
other i' ^ Jence angle (entry point ATM), and the brightness temperature translations 
(subrout > :(j TRANSREF), These major routines, and others called internally, are 
described and listed in Appendix A“1 . 


li 


SUBROUTINE ATM50 



EV- VERTICAL EMISSIVITY 
EH - HORIZONTAL 
EMISSIVITY 
AHEN - AHENUATION 
e- ACTUAL INCIDENCE 
ANGLE 

W - SURFACE WATER 
TEMPERATURE 
SAL - SALINITY 
TAPPV- VERTICAL 

APPARENT TEMPERATURE 
TAPPH- HORIZONTAL 
APPARENT TEMPERATURE 
DEFLT - DEFAULT VALUE 






ENTRY POINT ATM 


DEFLT 

ATTEN^ 



COMPUTE 

TRANSMIHANCE, 

AHENUATION 



EV - DEFLT 


CALLERS (VERT. POLARIZATION) 



EV- VERTICAL EMISSIVITY 
EH- HORIZONTAL EMISSIVITY 
AHEN- AHENUATION 
e- ACTUAL INCIDENCE ANGLE 
TW - SURFACE WATER TEMPERATURE 
TAPPV - VER. .CAL APPARENT TEMPERATURE 
TAPPH ■ HORIZONTAL APPARENT TEMPERATURE 
DEFLT - DEFAULT VALUE 









SUBROUTINE TRANS REF 



CALL DIECO/FRESNEL AT 
AT SURFACE WATER TEMPERATURE 
AND STANDARD INC IDENCE ANGLE 


EV- VERTICAL 
EMISSIVITY 
EH -HORIZONTAL 
EMISSIVITY 
e- ACTUAL INCIDENCE 
ANGLE 

STDANG - STANDARD 
INC IDENCE ANGLE 
SAL - SALINITY 
TW- SURFACE WATER 
TEMPERATURE 

VREF - VERTICAL REFERENCE 
BRIGHTNESS TEMPERATURE 
HREF - HORIZONTAL REFERENCE 
BRIGHTNESS TEMPERATURE 
DEFLT - DEFAULT VALUE 


HREF 

VREF 


DEFLT 




ll 
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APPENDIX A-1 


Description and LfsHng of Subroutines 

Fortran IV subroutines and a brief description of how to use them 
are listed in this appendix. Each routine is documented internally for more detailed 
user Information, The equations and formulas shown in the programs are based on 
derivations given in the text, with some modifications as to variable names. 

MAJOR SUBROUTINES 

Subroutine ATM50 A~1.I 

This routine computes the polarized emissivities and attenuation for the 
50° incidence angle case. Should input variables of incidence angle, water temp- 
erature, salinity, or vertical apparent temperature be default, the output 
variable will likewise be defaulted since all of these input quantities are necessary 
to compute emissivity and attenuation. Vertical emissivity and attenuation are 
computed independently of horizontal appc . ent temperature, so should horizontal 
apparent temperature be default, a default value only for horizontal emissivity 
will be returned along with computed values for vertical emissivity and attenuation. 

Entry point ATM A-1 ,2 

This entry point in subroutine ATM computes the polarized emissivities and 
attenuation for incidence angles other than 50°, The input and outj. ut arguments are 
identical with those in subroutine ATM50. The mean atmospheric temperature and 
the opacity, which are constant at all incidence angles, are determined in ATM50 
and are used here to determine the polarized emissivities and the attenuation. 

Subroutine TRANSREF A-1,3 

This routine performs the reference brightness temperature translations 
discussed in section II. Input quantities of water temperature, salinity, incidence 
angle, and standard incidence angle are required to calculate the Fresnel reflection 
coefficients and, likewise the emissivities, for a smooth sea. The values for 
frequency and salinity arc specified in a DATA statement. The change in emis“ 
sivity in translating from the surface water temperature at the actual incidence 
angle to the standard water temperature at the standard incidence angle is computed. 
Then, by knowing the emissivities for a roughened sea, the standard emissivities and 
the refe'^ence brightness temperature can be computed. 
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subroutine ATK5C ITH£T A , TW ,TA PV . TAPH .S AL . EV .EM, ATTENI 


s 

6 

7 

c 

c 

c 

8 

c 

9 

c 

10 

_c 


11 

12 

13 


THIS ROUTINE COMPUTES THE POLAPIZEO EMISSIVITIES ANO 
attenuation At 5t cegpees incident angle, the incident 
angle (TmETA) and water temperature <TW) are converted 
from degrees to radians ANO prom degrees centigrade 

TO degrees KELVIN INTERNALLY. THE ROUTINE ALSO CHECKS 
THE INPUT arguments FOR DEFAULT VALUES. SHOULD THERE 

EXIST any defaulted in«»ut parameters. The routine 

WILL RETURN CEFAJLT VALUES FOR THE CORRESPONDING 

outplt parameters.^ 


INPUT PARAMETERS 


1 " 
i «■' 

14 

C 

THETA = INCIDENT ANGLE (DEG» 

t-'i 

19 

C 

TM X WATER TEMPERATURE (DEG C) 


16 

C 

TAPV * vertical APPARENT TEMPERATURE (DEC K» 



i ■*. 

17 

C 

tafh X horizontal apparent temperature (OeG K» 

r j r* 

10 

C 

SAL X SALINITY IN PPM 

1 

19 

C 

output parameters 


20 

C 

ev * vertical emissivity 

I ' * 

21 

C 

EH X horizontal emissivity 


22 

C 

ATTEN X attenuation (0B» 

1 

23 

c 

INTEPNAL parameters 

1 ” 

24 

c 

1 X A FLAG hMICH is set to DETERMINE WHETHER OR NOT 

1 

25 

c 

ALL THE OUTPUT PARAMETERS ARE DEFAULTED FOR THE 

j 

26 

c 

50 DEGREE CASE, wmcn IsO, ALL THE 50 DEGREE CASE 

1 

27 

c 

PARAMETERS ARE OEFAULTEC ANC THERE IS NO NEED TO 

-• 

„2».^ 

. c 

PER»‘OKM ThP calculations AT ThE ENTRY POINT ATM. 


29 

c 

WHEN 1=1. ALl TmE OUTPUT PARAMETER VALUES FOR THE 


30 

c 

50 CEGREE CASE ARE NOT DEFAULT. I IS RESET TO 1 

X 

31 

c 

FOR each N-.W scan. 

i 

32 

c 

RV.RH X POLARIZED FRESNEL COEFFICIENTS 

7 : 

33 

c 

GAMMA s transmittance AT 5C DEGREES INCIDENT ANCLE 

- 

34 

c 

GmMANG X transmittance at OTHER INCIDENCE ANGLES 

M 

35 

c 

TA9ARV.TA3A3H s POLARIZED MEAN ATMOSPHERIC TEMP. 

)! 

36 

c 

TVAT.THaT = POLARIZED ATMOSPHERIC TEMPERATURES 

> 

37 

c 

T X ATMOSPHERIC OPACITY 

1 

38 

c 



39 


COMPLEX RV.RH 

•> 

40 


data :EFLT/03r7777777777/ 

- 

41 


DATA FREO. 1/13. 9.1/ 

1 

42 


IF (tmeta.eq.ceflt.or.tw.eq.deflt.or.tapv.eq.oeflt.or.sal.eo.deflt) 

! ■ . - 

43 


iGO TO 2C 


44 

c 

CONVERT incident AnGLE IN DEG 10 RADIANS 

- 

45 


ANGRA [xTHETA/37.3 

.1 

46 

c 

CONVERT WATER TEMPERATURE IN OEG C TO DEC K 


47 


TwK=TW^273.1« 

• 

48 

c 

CALCULATE THE VERTICAL EMISSIVITY FOR 

-* 

49 

c 

THE ZERO WIND CASE 


50 


CALL 0I£CC(FRE(J.TWK.SALI 

• 

51 


CALL FPESf.EL (ANGRAD.RV.RHI 

J 

52 


RtfM2=RV»C0NjG<PVI 
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53 EV«l.-«?VM2 _ _ 

54 C CALCULATE VERTICAL BRICHTNESS“TEHPERA TCRE I 

55 TBv«EV’TMK i 

56 C CALCULATE VERTICAL EXCESS TEMPERATURE 

57 TVEX=TAPV-T3V 

58 CALL £TTENV<TVEX, THETA, GAMHA,ATTEN,T> 

_ 59_ CALL TaTM (TVEX.CAHHA ,TM<,TA«>V,TVAT,TABARV» 

' 60 ' ‘ IF (TAPh.EQ.CtrLT) CO TO IP 

61 CALL HEXCESSITVEX.THEX) 

62 C l__.CALCULATE_MORIZONTAL BRIGHTNESS TEMPERATURE 

63 T3H*T4PH-t'hEX 

64 C calculate horizontal cMMiSIVlTY 

65_ Eh»TBH/THK 

66 ‘ CALL TATM J THEX .GAMMA .TWK'.fA PH~,fHAf,T A BARM » 

67 GO TO 30 

_ 68 IP EH«0EFLT 

69 CO TO 30 

7C 2C EV«0EFLT 

7l EHsOEFLT I 

72 ATTEN»0£FLT I 

73 I«0 

74 30 RETURN , 

75 ENTRY A TM (fnET A ,'th ,T a'pv , T APH ,EV, EH, A T TEN) 

76 C THIS ENTPy POINT CALCULATES THE POLARIZED ] 

77 _ C EMIS3IVITIES AND THE ATTENUATION FOR INCIDENT I 

78 C angles other THAN 5C DEGREES. THE MAIN ROUTINE 

79 C IS CALi.EC ONCE TS DETE'^^^INE THE POLARIZED MEAN | 

00 S _ atmospheric temperature <TA=ARV and TA9ARHI AnO THE ] 

01 C ATMOSP-ir.RIC opacity <T) for T“E 50 DEGREE CASE AND 

02 C • siNC- These arguments are constant at all other j 

83 C .incidence angles. THE EMISSIVITIES AND ATTENUATION ! 

84 C MAY D£ CETERMInED FUR ANY INCIDENT ANGLE. i 

85 C • • ! 

06 IFd.FQ.O. OS. THETA. EQ.DEFLTJ GO TO_70 i 

87 ANCRAD=ThETA/57. J . ! 

88 SECANG=1 ./COS(AnGRAD) \ 

89 c _ calculate the transmi tt ance_anq attenuation 

90 GAm'anC-=EXP«-SECANG*T» ' “ * t 

91 tTT£N=-i: .*410010 tOAMANG) 

92 IF (TAPV.EQ.DEFLT .OR.TH.EQ.CEFLTI, CO T0..40 ! 

93 THK=TW+273.18 ' 

94 call EPS (TAEARV.TAPV.GAMANG.THK.EVI ] 

95. . . GO TO 50 

96 4C EV«0EFLT ( 

97 5C IF(TAPh.EQ.0EFLT.3R.TH.£Q.DEFLT) GO TO 60 

96 , _ THK=TW*27’.18 

99 CALL EPS (TAPARH.TAPM.GAMANG.TWR.EHI 

103 GO TO 80 

_ 101 6C EH«DEFLT J 

102 GO TO 80 

103 70 EV»0FFLT 

104 Eh»DEFLT ^ 
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• 105 ATTEN»0EPLT 

r 106 ■ 00 RcTU«N ' 

I 107 EHO 

_ 1J^70 EQOtUlTY OR NOW-EQUfttlTY COMPAf<ISOM MAY NOT BE PE*NlNCFUL IN LOGICAL I^EXPRESSIOMS 
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• 1 

; i 

3 
_4 
» 
6 
7 
6 
9 

Ifi. 
11 
12 

13 

14 
19 
16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 
29 
3C 

31 

32 

33 

34 
39 

36 

37 
36 
39 
4C 

41 

42 

43 

44 

45 

46 

47 
46 
49 

90 

91 
52 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

Jl 

c 

c 

c 



c 

c 

c 

c 


c 


SUBROUT I TRftNSREFJTW, SAL .JMET4 , STDA NC, lEhtT^BV^Fj TBHP.EF »_ 

THIS ROUTINE CO'^RtCTS Trie RADIOMETER H£A3UREHENTS 

ev translating the brightness temperature to a 

reference 3R1GhT..£S3 TEMPERATURE BASED ON 293 DECREE 
KELVIN mater TeMPEf.ATURE AND THE STANDARD ANGLE 
£P«IFIED. 


INPUT FARAMZTERS 

_ TM«_ mater TiMPEPATURE_(OEG C> 

SAL • salinity in ppm 
TMETA s INCIDENT ANGLE COEG) 

STDAWG ■ STAN3ASD ANGLE JOEG> ! 

EV « vertical EMISSIVITY 
EH « horizontal emissivity 

OUJPUT PARAMETERS 

T8VREF « VERTICa“l REFERENCE TEMPERATURE TDEC K» 
TBMRiF « HORIZONTAL REFERENCE TEMPERATURE CDEC Kl 


DATA FREO.TS/13. 9.290./ . 

data OtFLT/0377777777777/ 

complex rw.rh 

IFTEV.rO.DEFLT.ANO. EH. EQ.DEF lTI GO TO 100* 

IF ITM.eQ.2EFLT.0R.SAL.eQ.DEFLT.0R. THETA. eQ.OEFLT. OR. STD ANG.EO. DEFL 
ITl GO TO ICO 

Convert Incident ancle in deg to radians 

AHGRADaTHETA/57.3 

CONVERT standard INCI CENT_ ANCLE 

IN DEGREES To RACIAt.S 


SAR»ST0AnG/57.3 

TMKsTM»273.l8 


CONVER T HA TER TEMPERATURE IN DEG C TO DEC K 

calculate the FRESNEL COEFFICIENTS ANO 
THE EMISSIVITIES AT Th£ CELL MATER 
temperature and TriE ACTUAL INCIDENT ANGLE 


CALL PIrCO(PPEQ. T;<<,EAL> 

CALL FFESNELCAn&paD.RV.RHI _ 

RVMZ«RV»CCNJG(RVI 

RHM2»KH*CCNJC(RH) 

EPSLV=1.-RVM2 

EPSLH«1. -RHM2 

calculate the Fresnel coefficients and 

the EMISSIVITIES AT THE STANOA^D WATER 

' ” temperature (29S. OeG KELVIN! AND THE 

STANCAKO incident ANGlE 

CALL DIECOtFPEQ. T3,SAL» ^ _ 

CALL FPESr.ELCSAR.RV.RH) 

RVMG2sRV*C0NJGCRV) 

RHMG2SP h*CONJG(RH) __ _ _ 

EFSLNV»1 . -PVMG2 
EPSLNn»l . -firiMG2 

calculate the Change in emissivity between 


70S74 Cl 06'g3-75 16t<»47 LASfL T 

53 C _ _ TM£ STAHQAH3 AUg C£U. WATER TEMPERATURE 

54 D£LEV»£PSLNV-EP’slV 

55 OtLEH»EPSLN«-EPSLM 

56 .C TK6NSLATE TO THE REFERENCE HATER TEHP. 

57 XFIEV.EQ.CEFLTI GO TO S3 

56 ESV>EV»OElEV 

55 TBVREF»TS*ESV • 

60 CO TO 60 . 

61 5C T3VREF«0EFLT 

67 60 IF (Eh.EO.DEFLTI G O TO SO 

63 ESH«Eh«DELEH 

64 T3HREP>TS«ESH 

^65 CO TO 7C0 

66 61 TBMREP*0EFLT 

67 CO TO 700 

66 100 T3VREF«0EFLT 

69 TBHBEF»DEFLT 

70 200 RETURA 





ROUTINES CALLED BY THE MAJOR ROLfTINES 


SubrouHne ATTENV A"1 ,4 

This routine which is called by ATM50 takes the 50” vertical excess tem“ 
perature c <d returns the attenuation and the transmittance* The incidence angl 
is then used in conjunction with the transmittance to obtain the atmospheric 
opacity. 

Subroutine HEXCESS A~1*5 

This routine, which is called at the entry point ATM, takes the 50® 
vertical excess temperature and returns the horizontal excess temperature. 
Subroutine TATM A"1 .6 

This routin , which is called at the entry point ATM, takes the input 
quantities excess temperature, transmittance, water temperature, and apparen. 
temperature, calc’. !atc<. the emissivity and atmospheric temperature and returns 
a value for the mean atmospheric temperature for any incidence angle. 
Subroutin e ** PS A"! .7 

This routine, which is c 'led at the entry point AIM, calculates the 
emissivity at any incidence angle given the mean atmospheric temperature, 
the apparent temperature, the transmittance, and the surface water temperature 
Subroutine FRESNEL A"1 .8 

This routine, called by both ATM50 and TRANSREK, computes the Fresnel 
reflection coefficients for sea water given the incidence angle. This routine 
must be initialized by calling entry point DIECO which computes the dielectric 
constant for sea water, given the frequency, water temperature, and salinity. 


I « 
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SUB»OUTINC ATTt’lV < TVLXSu , TmET* ,(.4?150 * ATT453 |T » 

TH!S KOUTINE CALCULATES THE ATTEHUATIOM, TR ANSNI 
ANO OPACITY at 5. CEGPCES ZNCIOCHT ANCLE. 

INPUT parameters 

TV£X53 « VERTiCAi. EXCESS TEMPERATURE 

_ AT 50 OiCREiS INCIOENT ANCLE lOEC 
TMETA » INCIOENT angle (DEC) 

OUTPUT parameters 

CAM50 » TRANSHITTANCi: AT 5C DEGREES _ 

ATTN30 ■ ATTENUATION AT 5b OtCRECS (OBI 
T ■ ATMOSPiERIC OPACITY 

X«TV£X50 

ATTn5:»:.:i795CS35*X-C.03005559*X**2*O.OC0004331*X**3 

.CAM5i»lC .••(-ATTN53/1C.I 

R/0»THrTA/57.3 

SEC>1./COS(PAOI 

T«-(AL0G(CAM5C)I/SEC 

RLTURI 

EHO 


TTANCC. 



I f 



<1 • 
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SOQROUTINE HExrCSS<V»Ml _ _ _ _ 

C* this routine CALCULiTEYtHEHORiZOMTAL EXCESS tYhPERATURE' 

C AS A function of vertical excess tekberature. 

C XMPUT'rtRAHEtER 

C V « VERTICAL EXCESS TENPERATURr COEC Kl 

C OUTPUT PARANPTER „ _ 

C H ■ horizontal EXCESS TEHPERATUKE (DEC X)" 

C 

H»l.53l32*V*.0(IALT»V**2-,a0CC»*V**3_ 

RETURK 

END 


LI 



I • 
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SUBROUT_INt TArn(T*X5C .gamma, TWK, TAPP, TATM50.TKI 

■ THISROUTINE calculates THE POLA^IZE3"ATMoSPMERIC' 
TEMPERATURE AT 50 DEGREES ANC MEAN ATMOSPHERIC TEMP, 


INPUT PARAMETERS 

TEX53 = EXCESS TEMPERATURE FOR 50 

TESREES INCIDENT ANGLE (DEC KI _ 

GAMMA s transmittance AT 5t DEGREES INCIDENT" ANGLE* 
THK 3 hater temperature (CEG K) 

TAPP = polarized apparent TEMPERATURE <OEC K) 


.1 

.j 

11 

12 

13 

C 

c 

c 

OUTPUT parameters 

TATM50 = POLARIZED ATMOSPHERIC TEMPERATURE (DEG K» 

TK = polarized mean atmospheric temperature 

i 

14 

c 



15 


E50=CTaPP-TEx50)/TMK 

i 

16 


A=TEX5C 

r 

1 

17 


b-gahma 

j 

10 


C = E50 


19 


d=thk 

f 

20 


TATM53=(A-(S-1.I *C*D-2. 7*B*9* ( 1 . -C ) ) / (3* (l.-C) ♦l.f 


21 


TK=TATM50/(1 ,-gammai 


22_ 


RETURN 


23 


END 

i • ^ 

5 


7 ] 
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subroutine EPSaK,TAPP,GAMANG.TMKtE» _ 

“This ROUTINE 'CALCULATES'THE POLARIZED EHISSIVI'TY AT A“ 
GIVEN INCIOENT ANGlE OThlR ThAN 50 DEGREES. 


1 

5 

V 

C 

INPUT PARAMETERS 

— — 


! 

6 

c 

t< = polarized mean atmospheric temperatlre 



, 

7 

c 

TAPP = POLARIZED APOARENT TtHPERATURr (DEG K) 



1 

9 

c 

GAMANG = transmittance 



j 

j 

9 

c 

TWK a water temperature <DEG K) 




1C 

c 

OUT?UT PARAMETER 




11 

c 

E a POLARIZED EMISSIVITY 




12 

c 





13 


TATM=TK*( 1. -GAMANG) 



1 

Ik 


A*TAPP 



1 

15 


B=GAMAN6 




16 


C=TATM 



1 

1 

ir 


DaTHK 



1 

18 


E= «A-C* (B*l.)-2.7*B*B)/Ca*(D-C-2.7*B) -C) 




19 


RETURN 




20 


END 




f- 


I I 
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• 

1 


SUBRCUTINP FRESMELCThiTi.RV.RM) 

350 

f. 

2 

C 

THIS routine computes ThE CLA*SSICAt 

3bC 

ij' 
1 1 

3 

c 

ERESfEL REELECTICN COEFFICIENTS FOR SEA WATER. 

370 


4 

c 

THIS SUGROUTINE MUST BE INITIALIZED BT 

38C 

« 

< 

5 

c 

CALLING THE SECONDARY ENTRY ROUTINE DIECO. 

390 

t 

6 

c 

THE complex PERMITTIVITY IS COMPUTED BY DIECO. 

4CC 

1 

7 

c 

ONCE HAVING computed THE PERMITTIVITY POR A PARTICULAR 

410 

'• 

8 

c 

FREOULNCV (FREQI , WATER TEMPERATURE < TEMP) , AND 

420 

1 

9 

c 

SALIN1TY(sAL> . FRESNEL MAY BF CALlEC REPEATEDLY 

43C 

1 

10 

c 

FOR DIFFERENT ANGLES (ThETAI . 

44C 


11 

c 


450 

" 

12 

c 

INPUT parameters 

46C 

i‘— 

13_ 

c 

FREQ = FcfQUENCY IN GHZ 

_470 

*: 

14 

c 

TEMP = TEMP£RATURE(DEG KJ 

~4 8C 


15 

c 

SAL = SALINITY IN PPM 

490 

'•1 

16 

c 

THETA s INCIDENT ANGLE (RAOIANSl 

SCO 

■f 

17 

c 

OUTPUT PARAMETERS 

510 


18 

c 

RV = VERT. POL. FRESNEL COEFF. (CCMPLFXI 

520 

•1 

19 

c 

RH = HORT. POL. FRESNEL COEFF, (COMPLEX) 

53C 

•1 

20 


OIMENSICI' E(2) 

54C 


21 


equivalence (Ed). ESP) 

550 

: 

22 


COMPLEX ESP. COST, SQ.RV.RH 

5tC 

-• 

23 


DATA f:D/3.031‘.15<)265/ 

57C 

•• 

24 


COST=CMPLX (COS (T 1ETA ), g.3 ) 

58G 

•* 

25 


SQ=CSORT (E3P-CmPlX(SIN ( ThET A ) **2 , 0 . 0 ) ) 

59C 

:t 

26 


RV»(ESP*COST-SQ)/(ESP»C03T*SQ) 

6C0 


27 


RH=(COST-SQ)/(COST»SO) 

610 

" 

28 


RETURN 

620 


29 


ENTRY OlECO(FP.EQ, TEMP, SAL) 

630 

7i 

30 

c 

■THIS Ef.TRY POINT CcMpUTcS THE DIELECTRIC CONSTANT 

640 

X 

31 

c 

POR SEA WATER. AlGO BASED ON PORTER’S WOPK (1R71). 

650 

- 

32 

c 

(SEE S.T. WU AND A.tC. FUNG'S REPORT NASA CR 2329). 

66C 

u 

33 

c 

. 

67'- 


34__ 

_ c_ _ 

CONVERT SALINITy.JO NORMALITY (SEE ST0GRYN_IEE£ _ 

680 

tf 

35 

c 

MTT AUG. 1971) 

69 c 


36 


S = SAL 

7CC 

> 

37 


Sa( (4.05 3E-09*S* 1,20 5E-35)*S^1. 707 e-02)*S 

- 7lC 


38 


T=TEMp-Z73.18 

720 

M 

39 

c 

normalized wavelength 

73 ; 


VO 


Rls((t.0C147*T-0.11)*T*3.3fi*(0.0l73*T-li.5 2) *S) • FREQ/30 .0 

740 


41 



75* 


42 


R2=(P1) **1.96 

76. 


43 


91=R1**.98 

77 


44 


E (2)=8T.C-15.3*S-0.363*T 

' 76 


45 


D»l.‘'*2.:*Rl*PI0»R2 

79 


V6 


£(1)=4.8+E(2)*(1.0*R1*PIO)/D 

60 


47 


SIG=(C.12'T.5.C)*S»C,04*T 

8: 


46 


E(2)=18. j*S:G/FREQ*E (2)*Rl/0 

82 


49 


RETUPN 

85 

•• 

50 


END 

86 


REPRODUCIBTT-TTY OF THE 
ORIGL\AJ. PAGE JS VOm 
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APPENDIX A-2 

Derivotion of the Mean Atm or oheric Temperature 

^SKY '* clefined an alrnospheric contribution to apparent temperature 
looking up through the atmosphere . ^SKY represented as the integral 

If the physical temperatures, . , 

e. 




T ( 2 ) is the physical temperature at altitude z the attenuation and 

the transmittance with "jr «» being the 

atmospheric opacity. Then, 

ix = 2 

Now, 1^5 i^Y ®><pressed as 

^To T ^ ec 0 


Define 




as 


,.To ^XiCc9 

( TCz)e. d 

Jo 


f Tq - X 5 ® 

Jo 

Using this value with the equation for T^j^y 


Integrating, 


TiXY = T^(G) [1 - 


Sample calculations from ATRAD show T^ to be virtually independent of 0 


30 
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